Casimir energy between a plane and a sphere in electromagnetic vacuum 



Paulo A. Maia Neto, 1 Astrid Lambrecht, 2 and Serge Reynaud 2 

'Institute) de Fisica, UFRJ, CP 68528, Rio de Janeiro, RJ, 21941-972, Brazil 
2 Laboratoire Kastler Brossel, CNRS, ENS, Universite Pierre et Marie Curie case 74, 
Campus Jussieu, F-75252 Paris Cedex 05, France 
(Dated: May 23, 2008) 

The Casimir energy is computed in the geometry of interest for the most precise experiments, 
a plane and a sphere in electromagnetic vacuum. The scattering formula is developed on adapted 
plane-waves and multipole basis, leading to an expression valid for arbitrary relative values of the 
sphere radius and inter-plate distance. In the limiting case of perfect reflection, the electromagnetic 
result is found to depart from the commonly used proximity-force approximation (PFA) significantly 
more rapidly than expected from scalar computations. 



The Casimir force is a mechanical effect of quantum 
vacuum fluctuations l| with a large impact in micro- and 
nanotechnology @, H | . Experimental advances allowing 
accurate measurements of small surface forces at micro- 
metric distances have led to a number of precise Casimir 
force measurements in the last 10 years (see reviews in 
[HHS])- These measurements and their comparison with 
Quantum ElcctroDynamics (QED) predictions have be- 
come a powerful tool for searching for the presence of 
the new hypothetical forces which are predicted by the 
models aiming at unifying gravity with quantum theory 
7, 8]. As the force varies substantially with the experi- 
mental conditions, accurate theoretical computations of 
realistic systems are needed for these comparisons to be 
reliable and fruitful. The influences of imperfect reflec- 
tion [|| and non null temperature (i, which have been 
studied extensively, are not discussed further here. 

Instead we will focus our attention on the rich connec- 
tion between Casimir effect and geometry As the 
most precise measurements performed to date involve a 
plane and a sphere, the effect of geometry is important 
for the purpose of theory-experiment comparison. It is 
usually calculated through the Proximity Force Approx- 
imation (PFA) [ll[ which amounts to average the plane- 
plane expression over the distribution of interplate dis- 
tances. This approximation can only be valid [13 , [l3| 
at the limit where the sphere R is much larger than the 
inter-plate separation L. Even in this limit, it is still 
worth specifying its accuracy in order to master the qual- 
ity of theory-experiment comparison [14| . 

The purpose of this work - A number of results going 
beyond the PFA have been obtained [IS EE El El El 
Sfl SI SI SI SI SE HI . It is only very recently that 
results were obtained for the case of direct relevance for 
the most precise experiments, namely the configuration 
of a plane and a sphere in electromagnetic vacuum [26| . It 
is the purpose of the present work to compute the Casimir 
energy in this configuration and draw consequences for 
theory-experiment comparison. 

We study the case of a sphere (radius R) on top of a 
plane plate (see Fig. 1). We denote C the center-to-plate 
distance and L the inter-plate distance. The plane-sphere 
Casimir energy £pg will be written in terms of a ratio to 
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FIG. 1: Plane-sphere geometry. 

the PFA formula (here written for perfect reflection) 
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Computations performed with scalar models [2J, |25| can 
be used to guess the expression of the ratio pps for small 
values of L/R. Considering that the electromagnetic case 
is given by the sum of results obtained with Dirichlet and 
Neumann boundary conditions, one thus gets 



PPS 



o 



L 2 
R 2 



(2) 



Bordag and Nikolaev [25[ remark that the electromag- 
netic result does not necessarily meet the scalar expres- 
sion but expect their magnitudes to be close to each 
other. The main result of the present work will be that 
the beyond-PFA corrections, as measured by the factor 
v, is in fact significantly larger in the electromagnetic 
case than in the scalar estimate. 

General scattering formula - Our starting point is the 
general scattering formula describing the Casimir energy 
between two scatterers in electromagnetic vacuum [6[. 
This formula has been used for dealing with rough or 
corrugated metallic mirrors in the limiting case where 
the roughness or corrugation was treated as a small per- 
turbation |1 SI SI- The PFA was thus found to be 
valid when the roughness or corrugation wavelength was 
larger than the other length scales. Here, we apply it to 
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the geometry of Fig. 1 where it is read 
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7^s a- n d 7^p are reflection operators on the spherical and 
plane scatterers while e represents one-way propa- 
gation along the z-axis between reference points sitting 
respectively at the centre of the sphere and on the plane. 
In the following, we write in a more explicit manner 
by introducing adapted mode basis. 

The plane-wave basis |k, 4>,p)^ (k the transverse 
wavevector, p = TE, TM the polarization, <fi = ±1 
the upwards/downwards propagation direction and £ the 
imaginary frequency) is well adapted to the description 
of free propagation and reflection on the plane: the oper- 
ator e~ lcc is diagonal in this basis (matrix elements e~ KC 
with k — y/£, 2 /c 2 + k 2 ); reflection on the plane also pre- 
serve all plane wave quantum numbers but <f> (the non 
zero elements of the matrix TZp are given by the Fresnel 
specular reflection amplitudes r p ). 

We also use the multipole basis \£mP)^ for the sub- 
space corresponding to a given £; £(£ + 1) and to are the 
usual angular momentum eigenvalues (with I = 1,2,..., 
to = —£,...,£) and P = E, M denote electric and mag- 
netic multipoles. By rotational symmetry around the 
z-axis, T> commutes with J z . Hence T> is block diagonal, 
and each block X)( m ) (corresponding to a given m) yields 
an independent contribution to the Casimir energy. The 
elements of this block £>( m ) are given by 
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When read from right to left, this expression has the fol- 
lowing interpretation: a multipole wave (£ 2 mP 2 ) is first 
decomposed into plane waves which propagate towards 
the (plane) plate where it is reflected; it then propagates 
back to the sphere, and is finally scattered into a new 
multipole wave (£\mP\). 

Reflection on the sphere can be written in terms of 
Mie coefficients ae(i£) and be(i£), corresponding to elec- 
tric and magnetic multipoles respectively [30| and of the 
finite rotation matrix [3l| elements S m i(9). The former 

depend on the reduced parameter £ = £R/c and the latter 
on the angle 9 such that cos 9 = ck/£ > 1, sin 9 = — «cfc/£, 
both evaluated for the imaginary frequency £. As soon 
as the Mie coefficients are computed, as the specular re- 
flection coefficients [32| , from the optical response of the 
sphere and plane, the scattering formula ((4]) allows one 
to obtain the Casimir energy for arbitrary separation dis- 
tances and arbitrary isotropic materials. 

Limit of perfect reflectors - In the present work, we fo- 
cus the attention on the case of perfectly-reflecting plane 
and sphere, which can be directly compared to already 



available scalar results. Precisely, we consider metallic 
materials described by the plasma model [12], with a 
plasma wavelength Ap. We then obtain a proper defini- 
tion of perfect reflectors as the limit where Ap is smaller 
than the length scales R and L. The Mie coefficients can 
thus be written in terms of the modified Bessel functions 
(defined as in (33j |) 

^ffi = Z[(-i)«+i <W(g-fe-i/ a (0 
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Meanwhile, specular reflection on the plane is described 
by ttm = — '"te = 1- In this case, we are able to sum 
up the TE and TM contributions in ((4| to get analytical 
expressions for the matrix elements 
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contain a factor exp(— 2££/c) which play a key role in 
the foregoing discussion of distance dependence of £ps • 

When the distance is larger than the radius, or equiv- 
alently when the sphere is small R <C C, the Mie co- 
efficients are needed for small 'size parameters' ( « 1. 
The^ dominant contribution thus comes from £ = 1 with 
oi (0 w -2&i(0 ~ -2£ 3 /3. The Casimir energy is then 
obtained from T>^ and T>^ , approximated by 2 x 2 ma- 
trices (P — E,M). The product of the nondiagonal ele- 
ments of these matrices is found to be negligible, and the 
following result obtained at lowest order in R/C 

9 fir R 3 

£ps= 16^1? ' [ A P< <i? «£] ( g ) 

Since the magnetic dipole contribution is one-half of 
the electric dipole one, the small-sphere limit ([5]) is 50% 
larger than the Rayleigh limit, calculated by accounting 
only for the electric contribution. Note that the latter 
can also be recovered from the general scattering for- 
mulas ([3]|3]), but not when using the Mie coefficients for 
the perfectly-reflecting sphere ([5]). As a matter of fact, 
perfect reflection corresponds to Ap being the smallest 
length scale whereas the Rayleigh limit corresponds to 
R <C Ap <C C, with a negligible magnetic contribution. 
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FIG. 2: Ratio pps showing the deviation from PFA as a 
function of L/R. The black and blue curves represent respec- 
tively the 'exact' result and the small-sphere limit (both for 
perfectly-reflecting plates) . 



As the sphere radius increases, higher values of £ 
and £ become increasingly important. The Mie coef- 
ficients grow as exp(2£) for £ ^> 1. When multiplied 
by the overlap integrals in ©, they produce a factor 
~ exp(— 2£L/c). Using the 'localization principle' [Hj], 
we may estimate the values of £ contributing appreciably 
to the Casimir energy for a given value of L/R. A given 
angular momentum £ corresponds semiclassically to an 
impact parameter B = c£/£. If B > R, its contribution is 
negligible since it corresponds to 'rays' that do not hit the 
sphere. With £ ~ c/L, we then expect that £ » R/L pro- 
vide negligible contributions. Numerical confirmations of 
this fact are presented below. 

Result of the numerical evaluation - We have numeri- 
cally evaluated the ratio pps of the plane-sphere Casimir 
energy normalized to the PFA expectation (see eqQJ 
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For perfectly-reflecting plates, this ratio is a function of 
the single parameter L/R, that we plot on Fig. 2. Sev- 
eral important features can be noticed on this plot. First, 
the small-sphere limit (blue curve) provides a good ap- 
proximation of the exact result (black curve) from large 
values of L/R down to L/R ~ 5. In the short distance 
limit L/R — > 0, the exact result goes to the PFA expec- 
tation pps = 1 whereas the small-sphere approximation 
leads to a strong underestimation. Note that we always 
have pps < 1, which means that the PFA systematically 
overestimates the Casimir energy. 

In order to discuss the most precise experiments, we 
zoom out on the interval L/R < 2 on Fig. 3, and de- 
vote a more detailed analysis to it. We first address 
the increased difficulty of the numerical evaluation near 
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FIG. 3: Zoom of Fig. 2 for L/R < 2. The green, blue and 
black curves are numerical evaluations of pps with maximum 
angular momentum t max = 10, 15, 27. The vertical dashed 
lines indicate the limit of validity of these evaluations. The 
red curve shows the scalar result (eq(2j at first order in L/R. 
The dots correspond to a best- fit quadratic function of L/R 
joining our numerical results to the PFA limit. 



the PFA (see [35| for a similar discussion in the geome- 
try of two equal spheres). As L/R is decreased towards 
the PFA limit, larger and larger values of I are needed. 
The localization principle requires I > a R/L where a 
is a numerical value that we have found to be approx- 
imately 4. The curves on Fig. 3 are numerical evalua- 
tions of pps respectively for maximum angular momen- 
tum £ max = 10,15, 27. The vertical dashed lines indicate 
the values L/R = 0.4 and L/R = 0.27 where the curves 
evaluated for £ max = 10 and 15 depart from the better 
calculation with ^ max = 27. Using the same argument, we 
predict that the curve computed with ^ max = 27 should 
be accurate down to L/R = 0.16 (also indicated by a 
dashed vertical line). It is worth doing computations 
with larger ^ max for increasing this range of validity [36| . 
But it also clear that these computations are expected 
to produce points close to the dots drawn by joining our 
numerical results to the PFA limit pps{L = 0) = 1. 

Comparison with existing results - In order to com- 
pare the present electromagnetic results with those ob- 
tained with scalar computations [24 , |25[ | , we have drawn 
the scalar prediction {2]), truncated at first order in 
L/R, as the red curve on Fig. 3. Clearly, it docs 
not fit the result we have obtained for the electromag- 
netic case. Precisely, the beyond-PFA correction appears 
to be several times larger in the electromagnetic case 
than in the scalar one. In order to make this remark 
more quantitative, we obtain a best-fit quadratic func- 
tion pps = 1 — v e xaL/R+ V2L 2 /R 2 joining the PFA limit 
pps(L — 0) = 1 to our numerical results. We thus get an 
estimate of the parameter v cra defined as in ^ but for 
the electromagnetic computation 
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This means that the Casimir energy between a plane 
and a sphere in electromagnetic vacuum departs from 
the PFA expectation significantly more rapidly than ex- 
pected from scalar computations. Our result agrees with 
that obtained in (26|, and it has been obtained through 
an independent calculation based on a different approach 
(i.e. the scattering approach of 0]). It constitutes a very 
important hint to be included in the discussion of the 
quality assessment of theory-experiment comparisons in 
the plane-sphere geometry |14| . 

This result has been obtained in the limiting case of 
perfect reflectors, for which some existing results were 
already available. Clearly, extra work is needed be- 
fore definitive consequences can be drawn for theory- 
experiment comparisons. As a matter of fact, the most 
precise experiments to date are performed at inter-plate 
distances not significantly larger than the plasma wave- 
length, and this entails that the effect of metallic response 



of the materials plays a non negligible role. However, the 
results already in our hands are a clear warning that the 
effect of plane-sphere geometry has to be treated with the 
greatest care when comparing the measured plane-sphere 
Casimir energy with QED theoretical predictions. 
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